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In the present paper we investigate the Tonks-Girardeau gas confined in a harmonic trap at finite 
temperature with thermal Bose-Fermi mapping method. The pair distribution, density distribution, 
reduced one-body density matrix, the occupations number of natural orbitals, and momentum distri¬ 
bution are evaluated. In the whole temperature regime the pair distribution and density distribution 
exhibit the same properties as those of polarized free Fermions because both of them depend on the 
modulus of wavefunction rather than wavefunction. While the reduced one-body density matrix, 
the natural orbital occupation, momentum distribution, which depend on wavefunction, of Tonks 
gas displays Bose properties different from polarized free Fermions at low temperature. At high 
temperature we can not distinguish Tonks gas from the polarized free Fermi gas by all properties 
qualitatively. 

PACS numbers: 05.30.Jp, 03.75.Hh, 67.85.-d 


I. INTRODUCTION 

For the high controllability and ’purity’ of quantum 
gas, it has been the popular platform to investigate novel 
quantum phenomena since it was realized in experiment. 
The cold atoms can be frozen to the zero point oscilla¬ 
tion in transverse with highly anisotropic magnetic trap 
or two dimensional optical lattices BHa, and thus the 
cold atoms become a one-dimensional (ID) quantum gas 
of strong correlation in longitudinal direction [3] . The in¬ 
teractions between atoms will be the effective ID contact 
interaction depending on the s-wave scattering length 
and its strength can be tuned from the weakly interact¬ 
ing regime to the strongly interacting regime with Fesh- 
bach resonance technique and confinement induced res¬ 
onance. In this way the Bose gas of strong repulsive 
interaction, Tonks-Girardeau (TG) gas [3 [3], can be re¬ 
alized, and offer researchers an opportunity to investigate 
the quantum many body physics of previously so-called 
’toy’ model. Besides the realization of arrays of ID gases, 
Jacqmin etc.’s experiment realized a single ID Bose gas 
close to the strongly interacting regime and the atom 
number fluctuation measurements can be performed in 
small slices of gases [7] . This allow us to study the prop¬ 
erties of ID Bose gas related with temperature and the 
interplay effect between interaction and temperature. 

Theoretically the unconfined ID Bose gas of contact 
interaction corresponds to the Lieb-Liniger model and 
can be exactly solved with Bethe-ansatz method in the 
whole interacting regime at zero temperature and at fi¬ 
nite temperature P P. Gombining the exact solution 
of integral model with local density approximation we 
can study the ID quantum gas confined in a harmonic 
trap BUHISl- The exact wave function of strongly inter- 
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acting TG gas in a harmonic trap can be obtained with 
Bose-Fermi mapping method [Sl ITilfTS] . So far most pre¬ 
vious research focus on the ground state properties of ID 
Bose gas. Its properties at finite temperature were not 
the focus of research although the temperature effect is 
important. The investigation of ground state properties 
show that as the atomic interaction becomes strong the 
ground state density profiles of Bose gas evolve from the 
Bose properties to the “fermionized” characters same as 
the polarized free Fermi gas [iniiiHiHi- At zero temper¬ 
ature, the equilibrium state of ID Bose gas is coherent 
quasicondensate in weakly interaction regime, and the 
gas reaches fermonized TG regime in strongly interact¬ 
ing regime. Although the TG gas show the same den¬ 
sity profiles as free Fermi gas, the density matrix and 
the momentum distribution exhibit the typical Bosonic 
properties mils]- 

Based on two-body correlation the interacting ID Bose 
gas can be classified as various physical regimes with the 
change of temperature and interaction and the 

interplay between interaction and temperature is an im¬ 
portant topic. At low temperature the ID Bose gas might 
be in the strongly interacting TG regime, a coherent 
regime or a fully decoherent quantum regime. Although 
the finite-temperature effect cannot be discounted, the 
theoretical determination of correlations, energy, and mo¬ 
mentum scales at finite temperature has only been par¬ 
tially investigated In the present paper we will 

investigate the TG gas trapped in a harmonic trap at 
finite-temperature. With thermal Bose-Fermi mapping 
method the pair distribution function, the reduced one- 
body density matrix (ROBDM), the occupation of natu¬ 
ral orbitals, and momentum distributions are evaluated 
for different temperature. As a comparison the corre¬ 
sponding result of polarized free Fermi gas are also ex¬ 
hibited. 

The paper is organized as follows. In Sec. II, we give 
a brief review of the model of ID TG gas and method. 
In Sec. Ill, we present the pair distribution function. 
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ROBDM, the occupation of natural orbitals, and momen¬ 
tum distributions of TG gases for different temperatures. 
A brief summary is given in Sec. IV. 

II. MODEL AND METHOD 


positive, zero, and negative, respectively. Since the wave- 
function of the polarized free Fermions is antisymmetry 
under exchanges, the constraint Eq. @ can be satisfied 
automatically. It is well known the many body wave- 
function of free Fermions can be formulated as a Slater 
determinants of single particle eigen function 4>u{x) 


In Bose quantum gases the properties of system are 
governed by the contact interaction between atoms, the 
strength of which depends on the s-wave scattering 
length as- When the cold atoms are confined inten¬ 
sively in transverse direction with tapping frequency uj ±_, 
the quantum gas becomes an effective one-dimensional 
many body quantum system. The effective ID interac¬ 
tion strength between atoms is 


giD 


4h'^al 

ma± 


(I-C 


Qs ^-1 
a± 


i^) 


with m being atom mass, the transversal length unit 
a_L = y/h/mui ±, and dimensionless constant C — 1.4603. 

For the ID quantum gas of N Bose atoms confined in 
a one-dimensional harmonic trap V{x) = with 

Lo being the trapping frequency, the system can be de¬ 
scribed by the Hamiltonian 


N 

^ = -E 
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_ 2m dx1 


V{xi) 


+ 9 lD X! ^{Xi-Xj), 


, ( 2 ) 

where Xi is the position of ith atom. Experimentally the 
effective ID interacting constant gio can be tuned by 
Feshbach resonance or confinement induced resonance in 
the full interacting regime from strongly attractive limit 
to strongly repulsive limit. 

In the present paper we shall focus on the strongly re¬ 
pulsive interaction limit, i.e., the TG gas. In this case 
the many body wavefunction ip^{xi,X 2 , ...jXn) is the so¬ 
lution of 


H'ip^{xi,X2 ,..., xn) = Ei;^{xi,X2, ■■■, Xn) (3) 


with H = - 


1 2 2 


For the strong 


2m dx^ 

interaction between atoms the wavefunction should sat¬ 
isfy the constraint 


i>^{xi, ...,Xi, ...,Xj, ...,xn) = 0 (4) 


a Xi = Xj {1 < i < j < N). By Bose-Fermi 
mapping method [HI ITS] to all the eigen wavefunction 
tj;^{xi,X2, ■■■,xn) of eigen equation Eq.([^ of Bose sys¬ 
tem there are corresponding wavefunction of polarized 
free fermions {x\^X 2 ^ ...jXat) 

lj^^{xi,X2, .■.,Xn) = A{Xi,X2, ...,Xn)iI^^{xi,X2, ■-,Xn), 

(5) 

where the function A is defined by A{xi,X 2 , ■■■,xn) = 
n sign(xi — Xj) and sign{x) is +1, 0, -1 for x being 


£ 1 

{xi,X 2 , ■■■,xn) = —i= det (j)i,i{xj). (6) 

vNi 

The eigen function of harmonic oscillator take the form of 
(l)u{x) = with H^{x) Hermite polynomial. 

As the set a = {vi, i' 2 , ■■■■, vn} take different combination, 
we can have all eigen wavefunctions including the ground 
state and excited states, which satisfies the eigen equa¬ 
tion a; 2 ,..., xat) = EN^ai’%^aixi,X 2 , ..^Xn), 

where the eigen energy Ejsf^a = + l/2)hoj. 

At zero temperature the system can be described by 
the ground state with a = {0,1,..., V — 1}, while at 
finite temperature the TG gases are governed by sta¬ 
tistical mechanical ensemble of Bosons. The TG gases 
are described by the probability PN,a attached to eigen¬ 
state i/jn aixi,X 2 ,..., a;Ar), which take the form of PN,a = 

N 

-vvith Z being the partition function, 

1=1 

P = l/ksT, the single particle energy = {vi + l/2)huj, 
and /i chemical potential. 

By the Bose-Fermi transformation Eq. © the eigen 
wavefunction of TG gas distinguish from the correspond¬ 
ing wavefunction of polarized free Fermions by the signa¬ 
ture function A(xi, X 2 ,..., xjv). Therefore the properties 
that depend on the modulus of wavefunction such as the 
pair distribution and single particle density distribution 
are same for TG gases and Fermi gases, and the proper¬ 
ties that depend on wavefunction such as ROBDM and 
momentum distribution are distinct with each other. 

Physically the pair distribution can be seen as the joint 
probability density that in two successive measurement 
if the first measurement finds one atom is at xi, then 
an immediate second measurement finds an atom at X 2 - 
At finite temperature the pair distribution function of 
TG gases and free Fermion gases are same, D 2 {xi,X 2 ) = 
D 2 {xi, X 2 ), which take the form of 


D 2 ixi,X2) = 


00 

E 

N=2,c 


P, 


N\ 


N,, 


{N-2)\J^ 


dxs... 


dxNi’M^ai’N^o 


Inserting the wavefunction into the above equation we 
can obtain the following simplified formula 

DI{xi,X2) = ^det^^piF {xi,xi) (7) 

2 

{1/1,1^2} P i=l 


where {pi,i' 2 } is the two dimension permutations of en¬ 
ergy level index {0,1,...,M} with M being the high¬ 
est single particle energy level to be considered. In the 
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calculation, we omit those high energy levels whose cor¬ 
responding Fermi-Dirac distribution fi, = is 

negligible small. 

The ROBDM can be viewed as the probability to 
find the particles at positions x and y in two successive 
measurements, respectively. At finite temperature the 
ROBDM p^{x,y) is defined as [T3] 


p\x,y) = 

N=1 a 


N! 

(N-n)! 



^2, •••, XN)'lpN,aiy’ ^ 2 , Xn)-{8) 


Here a sums over all possible permutations of N dimen¬ 
sion of {0,1,...,M}. In the calculations we also ignore 
the eigenstates with negligible fy. With Eq. (§ the 
Boson ROBDM p^{x,y) is expressed in term of the de¬ 
terminants of one-body Fermion density matrix pip^x, y) 
as follows [m |30] 


{x,y) 


Y — ^{sign{x-y)y dx2.-dxj+i 

^ det (9) 


With the properties of determinants dein,m P^ 
take the form of [30] 

det (10) 

l<n,m<j + l 

!^l...!^j + l P i=l 

Thus the multiple integral in Eq. can be reduced into 
single-variable integral, which makes the evaluation of 
ROBDM become possible. In the calculation more high 
single particle orbitals should be included in the sum¬ 
mation as temperature increases, so the computational 
amount of the exact implementation of Eq. ( [l0| will be¬ 
yond our ability. Eor high temperature we will follow the 
approximation of the one-body Eermion density matrix 
used in Eq. (10) in Ref. [50] . 

After we have reduced one-body density matrix p(x, y), 
which might be ROBDM of TG gas p^{x,y) or that of 
Fermions p^ {x, y), by solving the following eigen equation 
of p{x,y) 


p{x,y) (Px{y)dy = X^(pyx) 


( 11 ) 


we can obtain the natural orbital (pi{x) with Ai being 
the occupation number of natural orbital, which satisfies 
Ai = N. The momentum distribution can be evalu¬ 
ated by the Fourier transform of ROBDM 

/ OO 

dxdyp{x,y) (12) 

-OO 
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FIG. 1: (color online) Pair distribution for 5 Bosons and 
Fermions in a harmonic trap, which is same with each other, 
(a) T = O.lhuj- (b) T = 0.5hLo-, (c) T = l.Ohuj; (d) 
T = lO.Oftw. 


III. NUMERICAL RESULTS 


We will present the pair distributions, ROBDM, occu¬ 
pation number of natural orbitals and momentum distri¬ 
butions of ID TG gases in a harmonic trap with N = 5. 
As a comparison, we also show the results of polar¬ 
ized free Fermions. For simplicity in this paper we take 
and Tiu as the length unit and temperature unit. 

The pair distribution function D 2 {xi,X 2 ) depends on 
the absolute value of the wavefunction, so it is the same 
for the bosons and fermions. In Fig. [^the pair distribu¬ 
tion function of TG gas with TV = 5 are plotted at dif¬ 
ferent temperatures. At low temperature (for example, 
T = 0.1) the pair distribution function exhibit the same 
properties as those of TG gases at zero temperature and 
the temperature effect is not obvious [33]. The diagonal 
part (the region close to xi = X 2 ) means the probability 
that two Bosons appear at the same position, which is 
zero because the interactions between bosons are strong 
and two atoms separate as much as possible. With the in¬ 
crease of \xi — X 2 I, the joint probability density increases 
with oscillation first and decrease to zero at large enough 
distances. As temperature becomes high although the 
joint probability density also show the behavior that in¬ 
creases first and then decreases, the oscillation vanish 
completely. In addition, with the increase of tempera¬ 
ture the joint probability density D 2 {xi,X 2 ) decreases at 
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FIG. 2: (color online) Density profiles of the ground state of 
TG gases in a harmonic trap at different temperatures with 
A'’=5. The temperature is in nnit of Tiuj. 


small distance of \xi — X 2 I and distributes in larger re¬ 
gions. 

In Fig. I^the density profiles of TG gases with iV = 5 
at finite temperatures are plotted, which is the diagonal 
element of ROBDM p(x) = p^{x, x). At low temperature 
the density profiles display shell structure same as the 
case at zero temperature [T5] and Bose atoms occupy in 
the center region of harmonic trap. In this case atoms 
populate at lower single particle energy levels with large 
probabilities. As temperature become high the thermal 
energy of Bose atoms increases and the population prob¬ 
ability at higher energy levels will increase. The shell 
structure of density profiles vanish and Bose atoms dis¬ 
tribute in wider region. As temperature is high enough 
the property of Gauss distribution is exhibited. The 
density profiles can be fitted to the Gauss distribution 

The reduced one-body density matrix of TG gases 
p^{x,y) with Af = 5 at finite temperature are shown in 
Fig.i At low temperature the ROBDM display the sim¬ 
ilar properties to those at zero temperature. Although 
the ROBDM is diagonal dominant, the off-diagonal ele¬ 
ments are also not neglected. This is related to ODLRO 
and it is shown that there exists ODLRO in TG gases 
at low temperatures. With the increase of temperature 
p^{x,y) decreases more rapidly at large |a; — y\ than the 
diagonal part. At high temperature the off-diagonal ele¬ 
ments vanish completely and only the elements p’’{x,y) 
for small \x — y\ have finite values. As a comparison we 
display the ROBDM of polarized free Fermi gas p^{x,y) 
with = 5 in Fig. 1^ It is shown that ROBDM ex¬ 
hibit completely different behaviors from those of TG 
gases at low temperature. Even at very low temperature 
(T = O.lhui) the off-diagonal element of ROBDM are 
negligibly small and there exists not ODLRO. At high 



FIG. 3: (color online) Reduced one-body density matrix of 
TG gases in a harmonic trap at different temperatures with 
A=5. (a) T = O.lhuj-, (b) T = 0.5?ia;; (c) T = l.Ofttu; (d) 
T = lO.Ohoj. 


temperature the Fermi gases exhibit similar behaviors to 
TG gases. As the temperature is in the middle regime 
(for example, T = l.Ohuj) the off-diagonal elements of 
p^{x,y) are much larger thatn those of p^{x,y). 

The ODLRO is related to the occupation number of 
natural orbital Xi besides the behavior of ROBDM at 
large distance \x — y\. In macroscopic system if the occu¬ 
pation number Aq of the lowest natural orbital is much 
larger than one, the system will have ODLRO and exhibit 
BEC-like coherence effect. The system can be described 
by the associated natural orbital that play the role of 
an order parameter. In Eig. the occupation of nat¬ 
ural orbital of TG gases and polarized free fermions at 
different temperatures are displayed. Eor Fermi gas at 
T = Q.lhui the occupation number is one for the N low¬ 
est natural orbitals, while for higher natural orbitals the 
occupation numbers are zero. The distribution of occu¬ 
pations versus the orbital number i behave as a step-like 
function. As the temperature increase more higher natu¬ 
ral orbitals will be occupied and the occupation numbers 
of the lower natural orbitals decrease. The occupation 
distribution of TG gases at high temperature is similar 
to free fermions. But at low temperature the occupation 
of the lowest natural orbital is greatly larger than those 
of other naturals (Aq = 2.69 at T = O.lhuj). This means 
that the lowest natural orbital dominates the properties 
of the whole TG gases at low temperatures. With the in¬ 
crease of temperature the range of significantly occupied 
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FIG. 4: (color online) Reduced one-body density matrix of the 
ground state of polarized free fermions in a harmonic trap 
at different temperatures with N—5. (a) T = O.lhuj; (b) 

T = O.Shoj; (c) T = l.Ohuj; (d) T = lO.Oho;. 




FIG. 5: (color online) Occupation number of natural orbitals 
of TG gases (a) and polarized free Fermions (b) in a harmonic 
trap at different temperatures with N=5. The temperature 
is in unit of Uuj. 


higher orbitals increase and the lowest natural orbitals 
shall not dominate the properties of TG gases, we have 
to consider more natural orbitals to describe the system. 

In Fig. [§i the momentum distribution of TG gases 
with N = 5 are displayed at finite temperature. At low 
temperature momentum distribution behave as a sharp 
peak structure at zero momentum point, which is the 
typical property of Bose system. The Bose atoms appear 
around the zero momentum region as the largest prob¬ 
ability and the occupation probability decrease to zero 
rapidly at finite momentum regime. The height of peak 


FIG. 6: (color online) Momentum distribution of TG gases 
(a) and polarized free Fermions (b) in a harmonic trap at 
different temperatures with N=5. The temperature is in unit 
of hu!. 


shrinks with the increase of temperature and Bose atoms 
occupy large momentum regime because of the assistance 
of thermal energy resulting from the finite temperature. 
At high temperature the sharp peak structure of momen¬ 
tum distribution disappears and Bose atoms distribute in 
wider momentum regimes. As a comparison the momen¬ 
tum distributions of polarized free Fermions of = 5 
are displayed in Fig. [6)3. The momentum distributions in 
momentum space are exactly same as the density profiles 
in real space. It deserves to emphasize that the difference 
from those of TG gases mainly exhibit at low tempera¬ 
ture. In this situation free Fermi gases show shell struc¬ 
ture same as the distribution in the real space. This is 
the main difference between the Bose system and Fermi 
system. Similar to the properties of ROBDM we also can 
not distinguish TG gases and polarized free Fermi gases 
by the momentum distributions at high temperature. 


IV. CONCLUSION 

In conclusion with thermal Bose-Fermi mapping 
method we calculated pair distribution function, density 
profiles, ROBDM, the occupation number of natural or¬ 
bitals, and momentum distribution of strongly interact¬ 
ing TG gases at finite temperature. 

In the whole temperature regime the diagonal part of 
pair distribution function D 2 {xi,X 2 ) is zero, which shows 
that two strongly interacting atoms always separate as 
much as possible. Gomparison to the case at low tem¬ 
perature D 2 {xi,X 2 ) distribute in larger regimes at high 
temperature . At low temperature the density profiles 
exhibit the Fermion-like shell structure, and the shell 
structure vanish with the increase of temperature. At 
high temperature the density profiles can be fitted as 
Gaussian distribution p{x) = ^=^exp(—i^). The in¬ 
vestigation of ROBDM shows that at low temperature 
there exist ODLRO in TG gas similar to the properties 
of zero temperature case, while the ODLRO vanish at 
high temperature. Gorrespondingly, at low temperature 
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momentum distribution displays a sharp peak structure 
as if many Bose atoms ’condensate’ at zero momentum. 
While with the increase of temperature the sharp peak 
disappear and TG gas distribute in high momentum re¬ 
gions. 

The pair distribution and density profiles depend on 
the modulus of eigenfunction so they are same for TG 
gas and polarized free Fermions in the whole tempera¬ 
ture region. The ROBDM and therefore the occupation 
number of natural orbital and momentum distribution 
depends on the eigenfunction so their properties for TG 
gas are different from those of free Fermions. But the 


distinction are shown only at low temperature, and we 
can not distinct them from each other qualitatively at 
high temperature. 
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